Theory of smeared quantum phase transitions 
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We present an analytical strong-disorder renormalization group theory of the quantum phase 
transition in the dissipative random transverse-field Ising chain. For Ohmic dissipation, we solve 
the renormalization flow equations analytically, yielding asymptotically exact results for the low- 
temperature properties of the system. We find that the interplay between quantum fluctuations and 
Ohmic dissipation destroys the quantum critical point by smearing. We also determine the phase 
diagram and the behavior of observables in the vicinity of the smeared quantum phase transition. 
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One of the most basic questions concerning phase tran- 
sitions in random systems is whether or not a sharp tran- 
sition survives in the presence of quenched disorder. Ini- 
tially, it was suspected that disorder destroys any critical 
point because different spatial regions order at different 
temperatures. However, it was soon realized that classi- 
cal continuous phase transitions generically remain sharp 
in the presence of weak disorder because finite spatial re- 
gions cannot undergo a true phase transition (see Ref. [1] 
and references therein). 

Nonetheless, rare strongly coupled spatial regions play 
an important role. They can be locally in the or- 
dered phase even if the bulk system is in the disordered 
phase. The slow fluctuations of these regions give rise 
to a singular free energy in a whole temperature re- 
gion around the transition (called the Griffiths phase) 
[2, 3]. In generic classical systems, this is a weak ef- 
fect, because the Griffiths singularity is only an essen- 
tial one. In contrast, rare regions can play a more im- 
portant role at zero-temperature quantum phase tran- 
sitions where order-parameter fluctuations in space and 
(imaginary) time need to be considered. Quenched dis- 
order is perfectly correlated in time direction, and this 
enhances the Griffiths singularities. In the prototypical 
random transverse-field Ising systems, the singularities 
take power-law forms, implying, e.g., a divergent suscep- 
tibility in the Griffiths phase [4, 5, 6]. The transition 
itself is governed by an exotic infinite-randomness criti- 
cal point [7, 8], but remains sharp. 

Recently, it was noted that dissipation can further en- 
hance rare region effects at quantum phase transitions 
with Ising order parameter symmetry. Each locally or- 
dered region acts as two-level system. When coupled 
to an Ohmic dissipative bath, it can undergo the local- 
ization transition of the spin-boson problem [9]. Thus, 
each region can order independently of the bulk system, 
destroying the sharp phase transition by smearing [10]. 
In view of this observation, it would be highly desirable 
to treat the nonperturbative physics of these dissipative 
rare regions within the framework of the renormalization 
group (RG) commonly used to describe phase transitions. 
Such a theory would not only unveil the ultimate fate of 
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Figure 1: (Color online) Zero temperature phase diagram and 
magnetization of the dissipative random transverse-field Ising 
chain as a function of the typical transverse field ht yp . SO 
and SD denote the strongly ordered and disordered conven- 
tional phases; WO and WD are the weakly ordered and dis- 
ordered quantum Griffiths phases, (a) No dissipation: sharp 
QCP. (b) Ohmic dissipation: smeared transition with the in- 
homogeneously ordered (IO) phase replacing the WD Griffiths 
phase, (c) Distributions of the bonds J and fields h e s in the 
various phases. The shaded area quantifies the fraction w of 
J's bigger than h e s's [see Eq. (10)]. 



the critical point, it would also predict quantitatively the 
behavior of many observables near the transition. 

An important step towards this goal was taken by 
Schehr and Rieger [11, 12] who studied the dissipa- 
tive random transverse-field Ising chain by a numerical 
strong-disorder RG. They confirmed the smeared tran- 
sition scenario and focused on the infinite-randomness 
"pseudo'-critical point arising at intermediate energy 
scales where dissipative effects are less important. 

In this Letter, we develop a comprehensive strong- 
disorder RG for the dissipative random transverse-field 
Ising chain. We derive RG flow equations for the distri- 
butions of the fields, bonds and magnetic moments and 
solve them analytically, providing asymptotically exact 
low-energy results. We prove that the quantum critical 
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point (QCP) is destroyed by Ohmic dissipation. Instead, 
a smeared quantum phase transition separates a conven- 
tional paramagnet from an inhomogeneously ordered fer- 
romagnet (Fig. 1). In the remainder of the Letter, we 
sketch the derivation of our theory, compute important 
observables, and discuss the relevance of our results. Ex- 
tensive details will be given in a longer paper. 

Our starting point is the dissipative random 
transverse-field Ising chain defined by the Hamiltonian 

i i 
i,n i,n 

where a*' z are Pauli matrices. The bonds Ji and fields 
hi are independent random variables; a\ n (a^ n ) are the 
creation (annihilation) operators of the n-th oscillator 
coupled to spin Oi via \i, n , and v^ n is its frequency. Ini- 
tially, all baths have the same Ohmic spectral function 
£((J) = nY, n A L^( W ~ v i,n) = 27rawe-"/"<=, with a the 
dimensionless dissipation strength and u c the (bare) cut- 
off energy. (The cutoff will change under the RG and the 
dissipation strength will become site-dependent.) 

To characterize the low-energy behavior of the system 
(1), we now develop a strong-disorder RG [13, 14]. The 
idea of this method is to successively integrate out local 
high-energy modes. In our case, the competing energies 
are the transverse fields, bonds, and oscillator frequen- 
cies. Each RG step proceeds as follows: We first find 
the largest energy in the system ft = max(/ii, J,;, uj c /p) 
where p S> 1 is an arbitrary constant [15]. We then lower 
the energy scale from f2 to fl — dfl by (i) integrating 
out all oscillators (at all sites i) with frequencies between 
p(fl — dO) andpO and (ii) decimating all transverse fields 
and bonds between (Q — dfl) and ft. 

For j; > 1, the oscillators can be treated using adi- 
abatic renormalization [9]. As a result, the transverse 
fields renormalize according to 



K = hi cxp 



p(n-dfi) u 



while the bonds remain unchanged. Here ai is the renor- 
malized dissipation strength at site i. 

To decimate a strong bond Ji = fl, we assume the spins 
Oi and <Ji+i to be locked together as an effective spin 
cluster a with moment /x and renormalized transverse 
field h obtained in second order perturbation theory, 



h = hih i+x /Ji . 



(3) 
(4) 



a couples to a renormalized bath of dissipation strength 

a = a t + a l+ i = a(jii + fi i+ i) = ap, . (5) 



For a strong field, hi — fi, the corresponding spin Oi 
is delocalized in a z basis and thus eliminated, creating a 
new bond between sites i — 1 and 



J — Ji-iJi/hi 



(6) 



Note that for spins about to be decimated, hi = 
Q is the fully renormalized tunnel splitting hi = 
hioiphiQ/ujco) 01 ^^ 1 ^ ^ where hio and lo c q are the field 
and bath cutoff of the i— th cluster when it was formed 
at the higher energy uj c q/p- 

The recursion relations (3), (4) and (6) are identical 
to the dissipationless case [7], the baths enter only via 
(2) together with the renormalization of the dissipation 
strengths (5). Our RG procedure is related to the one im- 
plemented numerically by Schehr and Rieger [11]. How- 
ever, treating the oscillator modes on equal footing with 
the other degrees of freedom (by reducing the bath cut- 
off globally in each step) allows us to solve the problem 
analytically. 

The complete RG step consisting of recursion relations 
(2)-(6) is now iterated with the energy scale Q being de- 
creased. At each stage, the remaining bonds J and fields 
h are independent, but the fields and magnetic moments 
are correlated. Using logarithmic variables T = ln(fij/fi) 
[where fli is the initial (bare) value of fi], C = In(f2/J) 
and (3 = ln(f2//i), we can thus derive RG flow equations 
for the bond distribution 'P(C) and the joint distribution 
of fields and moments lZ(f3, y). They read 



dV 
df 

dK 
9f 



&P ( C 

— + (l-aPo)ft/j(0) \V®V 

+ [V(0)-(l-ap )K f3 (0))V , (7) 

(i-a^) — + v(o) [n ® nj 

-{V(0)-(l- a jI Q )K p (0)}K , (8) 



where 1Zp{(3) = / °° 72.(/3, A*)d/it is the distribution of the 
fields and ~p is the average moment of clusters about to 
be decimated (defined by ~p 7lp(0) = Jl°° fiR(0, ju)d/x). 

The symbol V g> V = j^V(C')V(C - C') d C denotes the 
convolution. The first term on the r.h.s. of (7) and (8) 
is due to the rescaling of C and f3 with T and the renor- 
malization (2) of h by the baths. The second term im- 
plements the recursion relations (3), (4) and (6) for the 
moments, fields and bonds. The last term ensures the 
normalization of V and 1Z. As expected, for a = 0, (7) 
and (8) become identical to the dissipationless case [7, 8]. 

Important insight can already be obtained from the 
structure of the flow equations. The probability of deci- 
mating a field, (1 — a]I )7lp(0), decreases with increasing 
dissipation strength and cluster size. Clusters with mo- 
ment jU > 1/a are not decimated. Thus, in the presence 
of dissipation, the flow equations always contain a finite 
length scale above which the cluster dynamics freezes. 
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We now search for stationary solutions of the flow 
equations (7) and (8) that describe stable phases or crit- 
ical points. There are two trivial cases: If all bonds are 
larger than all fields, only bonds are decimated, build- 
ing larger and larger clusters. This is the conventional 
strongly ordered (SO) ferromagnetic phase. If only fields 
are decimated, we are in the conventional strongly disor- 
dered (SD) paramagnetic phase. 

In the more interesting case of overlapping field and 
bond distributions, we look for solutions invariant under 
a general rescaling 77 = (/f((T), 9 = (5/ fp(T) and v = 

Without dissipation, a = 0, there are three types of 
well-behaved solutions [8] : a line of fixed points (parame- 
terized by 7^o) with fp = l,f(= cxp(lZoT) and the aver- 
age moment increasing as T. It corresponds to the weakly 
disordered (WD) Griffiths phase. There is another line 
of fixed points with = 1 and fp = / M = cxp(Vo^) 
(parameterized by Vq) which corresponds to the weakly 
ordered (WO) Griffiths phase; and, separating these two 
phases, an infinite-randomness QCP with fa = fp = T 
and fn = with 2(f> = 1 + a/5. 

In the presence of dissipation, a ^ 0, the scenario 
changes dramatically. For overlapping bond and field 
distributions, we found only one line of well-behaved 
fixed points (parameterized by Vq > 0) corresponding 
to the ordered phase [16]. Here, fa = 1, = exp('Por), 
fp = rcxp('Por). The fields become much smaller than 
the bonds, justifying the perturbative treatment of the 
RG step. The fixed-point distributions are 

V*(0 = Voe~ v ^ , (9a) 
n*{6,v) = Uoe-^Sie - av) , (9b) 

i.e., fields and moments are perfectly correlated. Here, 
IZo is a nonuniversal constant. This fixed point is similar 
to the WO Griffiths phase for a = 0, but fp/ — > 00 as 
r — > 00. Transforming the field distribution (9b) back to 
the original transverse fields h gives power-law behavior 
~ h 710 ^"^^ 1 . We could not analytically solve for the 
nonuniversal constants Vo and TZq in terms of the bare 
distributions and a. Their numerical values will be given 
elsewhere. 

We emphasize that we have shown that there is no 
fixed point solution with /<; / fp — > const as V — > 00 in the 
presence of dissipation, implying that there is no QCP 
where fields and bonds compete at all energy scales. This 
important result proves that Ohmic dissipation destroys 
Fisher's [7, 8] infinite-randomness critical point. Physi- 
cally, it is due to the fact that finite spin clusters (of size 
~ 1/a) can develop true magnetic order. 

The complete low-energy thermodynamics can be ob- 
tained from the RG fixed point solutions. To characterize 
the phase diagram (Fig. 1) in terms of the bare variables 



we introduce the probability 

roc pj 

w= dJPi(J) dh eS Rj(h eS ) , (10) 
Jo Jo 

of a bare bond J being greater than an effective field 
(a bare field, fully renormalized by the baths) h e g = 
h{ph/u c ) a / { ' 1 - a \ P/(J) and Ri{h eS ) are the bare initial 
distributions of these variables [see Fig. 1(c)]. 

For w = and w = 1, these distributions do not over- 
lap. The system is in one of the conventional phases (SD 
or SO) without Griffiths singularities where disorder is 
RG irrelevant. For < w <C 1, arbitrarily large rare 
clusters can form under renormalization. Without dissi- 
pation, a = 0, these clusters have small but nonzero effec- 
tive fields. They thus slowly fluctuate, and the system is 
in the WD Griffiths phase [see Fig. 1(a)]. In the presence 
of dissipation, a/0, clusters with moment (X > 1/a have 
zero effective field. They freeze and order independently 
from the bulk. The sharp transition is thus destroyed by 
smearing, and the WD Griffiths phase is replaced by an 
inhomogeneously ordered (10) ferromagnetic phase [see 
Fig. 1(b)]. Finally, with w approaching 1, the system de- 
velops bulk magnetic order but rare fluctuating clusters 
still exist, i.e., we are in the WO Griffiths phase. In the 
presence of dissipation, the IO and WO phases are sepa- 
rated by a crossover rather than a QCP. The asymptotic 
low-energy properties of both phases are described by the 
solution (9) with Vo monotonically decreasing with w. 

We now turn our attention to observables near the 
smeared phase transition, focusing on the IO ferromag- 
netic phase which is the novel feature of our system. The 
magnetization is dominated by the large frozen droplets 
which arise in rare regions where the bonds are greater 
than the local fields. Because they are static, any weak 
coupling mediated by the bulk is sufficient to align them. 
Hence, the magnetization is proportional to the volume 
of the rare frozen droplets, which for a and w <C 1, is 

TO - W 1/a . (11) 

The low-temperature magnetic susceptibility can be 
computed by running the RG to energy scale fl = T 
and assuming the remaining spin clusters to be free. For 
asymptotically low energies, the RG flow is dictated by 
the fixed-point solution (9), leading to 

X-T- 1 - 1 ^, (12) 

with z = \ jVo- Note, however, that at higher energies, 
the flow is dominated by strong fields and the suscepti- 
bility therefore behaves as in the weakly disordered un- 
damped Griffiths phase [8]: 

X - <5 4 ~ 20 [In (1/T)] 2 T- 1+1 / z ' , (13) 

with z' w 1/(2(5), and 5 ~ (ln/i e ff) — (In J). The crossover 
energy fl c separating the two regimes can be estimated 
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as the energy in which the high-energy mean moment 
cluster, ~JL ~ TS 1- ^, reaches the critical size 1/a. Hence, 
aln(0//f2 c ) ~ t^ -1 . Below fi c , the mean magnetic mo- 
ment increases much more rapidly, ~p ~ exp (VoT). 

In summary, we have developed an asymptotically ex- 
act strong-disorder RG theory for the dissipative random 
transverse-field Ising chain. We have solved the result- 
ing flow equations analytically and proven that the QCP 
is destroyed by smearing. The smearing is the result of 
the interplay between disorder and dissipation. Dissipa- 
tion alone leads to a conventional critical point [17], while 
disorder alone leads to an exotic infinite-randomness crit- 
ical point [7, 8], but the transition remains sharp. In the 
remaining paragraphs, we put our results in broader per- 
spective, and we discuss further implications. 

We first consider the dissipative random transverse- 
field Ising model in higher dimensions. The recursion 
relations (2)-(5) are the same as in one dimension while 
decimating a field now generates couplings between all 
nearest neighbor sites, changing the topology of the lat- 
tice. An analytical solution of the RG flow equations 
thus appears impossible. However, the dissipation terms 
are local and take the same form as in one dimen- 
sion. In particular, the probability of decimating a field, 
(1 — ap )TZ/3(0) , is reduced with increasing dissipation 
and vanishes for clusters with finite moment /U > 1/a. 
Thus, a critical fixed point solution is impossible, and 
the infinite randomness critical point found in the dissi- 
pationless case [18, 19] is destroyed by smearing. More- 
over, the weakly disordered Griffiths phase is replaced 
by the inhomogeneously ordered ferromagnet. Note that 
Ohmic dissipation also suppresses the quantum Griffiths 
singularities at the percolation quantum phase transition 
[20] in a diluted transverse-field Ising model [21]. How- 
ever, the percolation transition remains sharp because it 
is driven by the critical geometry of the lattice. 

Our results for a dissipative Ising magnet must be 
contrasted with the behavior of systems with continuous 
O(N) symmetry. While large Ising clusters freeze in the 
presence of Ohmic dissipation, 0(7V) clusters continue 
to fluctuate with a rate exponentially small in their mo- 
ment [22]. This leads to a sharp transition controlled by 
an infinite-randomness critical point in the same univer- 
sality class as the dissipationless random transverse-field 
Ising model [23]. All these results are in agreement with 
a classification of weakly disordered phase transitions ac- 
cording to the effective dimensionality of the rare regions 
[24]. If their dimension is below the lower critical dimen- 
sion d~ of the problem, the behavior is conventional; if 
it is right at d~ , the transition is of infinite-randomness 
type; and if it is above d~ , finite clusters can order inde- 
pendently leading to a smeared transition. 

To the best of our knowledge, this work is the first 
quantitative analytical theory of a smeared phase transi- 



tion. The results directly apply to quantum phase tran- 
sitions in disordered systems with discrete order parame- 
ter symmetry and Ohmic damping. Our renormalization 
group approach should be broadly applicable to a variety 
of disordered dissipative quantum systems such as arrays 
of resistively shunted Josephson junctions [25, 26]. 
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